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• Abstract 

Explicit counterterm construction is presented in a symmetry breaking background for the 
2PI effective action of the self-interacting one-component scalar field including skeleton dia- 
grams to C(A 2 ). The applied strategy is an extension of our treatment of the 2PI-Hartree 
approximation [24] . The procedure is also applied to the O(N) model. 
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Renormalisation invariance is a powerful concept in constructing field theoretical models describing 
phenomena in a restricted range of momenta. It imposes insensitivity requirement for the physical 



0© predictions to the choice of the maximal momentum if the range characterising the phenomenon 
^ ' under investigation is much below that value. Renormalisable models represent a specific subclass 
where the maximal momentum could be sent to infinity. In view of triviality of scalar field theories, 
^ ■ renormalisability just ensures that high precision predictions calculated with finite (but large) cut-off 
in such models will not depend sensitively on the cut-off. 

Renormalisability criteria for perturbatively computed physical quantities as power series of the 
renormalised couplings represent a classical piece of knowledge [T]. One constructs at every order 
also the corresponding approximation of the counterterm Lagrangian, which cancels the cut-off de- 
pendence of some loop integrals. The renormalisation of partially resummed perturbative series is 
nontrivial because in association with the selected subset of the perturbative diagrams one has to 
perform also the resummation of appropriately chosen counterterm diagrams. 

In the past 6-7 years considerable progress has been achieved in constructing renormalisable two- 
particle-irreducible (2PI) approximate solutions of perturbatively renormalisable theories. The 2PI 
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approximation was developed into a valuable non-perturbative tool in studies of collective features 
of quantum field theories both in and out of equilibrium. With its help one can investigate the 
evolution of the spatial distribution of ground state field condensates as well as propagators of their 
quasi-particle excitations. 

In a series of papers [21 El E] van Hees and Knoll proved renormalisability of the 2PI approximate 
solutions in the real time formulation of quantum field theories with help of the subtraction proce- 
dure designed by Bogoliubov, Parasiuk, Hepp, and Zimmermann (BPHZ). Their study focused on 
the symmetric phase with no condensate in the ground state. This line was continued in 2003-04 
by Blaizot, Iancu and Reinosa [51 [6] having demonstrated iteratively that self-consistently solved 
propagator equations resum infinite subsets of different orders of the perturbation theory. Mainly 
based on diagrammatic analysis of the iterations they also constructed the counterterms necessary 
to balance the cut-off dependence of the loop integrals appearing in the self-consistent equations. 
Moreover they pointed out the importance of the large momentum asymptotics of the propagator for 
the determination of the cut-off dependence of the counter couplings. Related investigations were 
done also using the auxiliary field formulation of <£> 4 theories (e.g. after Hubbard- Stratonovich trans- 
formation) [3, El E] and also concentrating on the problem of gauge fixing parameter dependence of 
2PI approximation in gauge theories (TQl HI] • 

A complete description of how to renormalise arbitrary n-point functions derived in 2PI-formalism 
was presented by Berges, Borsanyi, Reinosa and Serreau in Ref. [12], also for non- vanishing expecta- 
tion value of the scalar field. The detailed analysis of that paper was mainly based also on diagram- 
matics and explicit counterterm formulae were displayed only partially. Since then several numerical 
implementations were reported for scalar theories [T3l HU [15] and also for QED [T6l Ej- The sym- 
metry aspects of the 2PI- formalism was investigated in [4j. It was shown that only the so-called 
"external N-point functions" , i. e. derivatives of the effective action with respect to the mean field 
fulfill Ward identities. In particular, the self-consistent propagator differs from the external two-point 
function and in consequence it violates Goldstone's theorem. The validity of Ward identities was 
discussed recently in the context of QED in |18j . 

Investigations of the renormalisation were performed also for a related resummation scheme us- 
ing local (momentum-independent) insertions, the so-called two-particle-point-irreducible (2PPI) ap- 
proximations. The renormalisation of this scheme was proved using counterterms [19] and the scheme 
was used in numerical works both in equilibrium [20] and out of equilibrium [2U [22]- Recently in 
Ref. |23j a renormalisation group inspired approach was developed to investigate the connection be- 
tween momentum-dependent resummation and renormalisation. There a comparison of the method 
with renormalisation methods relying on Bethe-Salpeter-like equations was also given. 

The purpose of the present paper is to deduce the resummed counterterms in a transparent way, 
and provide explicit formulae as far as possible for the counter couplings. In a previous paper [21] 
we have shown for the 2PI-Hartree approximation that results of the iterative renormalisation of 
Ref. [6] can be obtained also with help of a one-step substitution. In this paper we demonstrate 
first the applicability of our procedure to all cases where the large momentum asymptotics of the 
propagator remains the same as it was in the tree level approximation. In particular, we find that in 
such cases there is no need to use any auxiliary Bethe-Salpeter-like equation for the determination 
of the 4-point counter-couplings. This feature can be studied neatly on the example of the complete 
two- loop skeleton diagram set of the 2PI effective action in a constant, homogeneous symmetry 
breaking background. The basic definitions and notations of the two models considered in this paper 
are introduced in Section 2. In Section 3 the construction is performed first for the case of the 
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one-component real field. It is applied to the more involved case of the O(N) model in Section 5. 
The presentation will be detailed and still more compact than was the case of previous analyses. 

The generalisation of our procedure to the case of the modified large momentum behaviour of 
the propagator is discussed for the one-component scalar field in Section 4. It leads to a new 
important intermediate step in designing the algorithm of the counterterm determination. Namely, 
the separation of the divergent pieces of different skeleton diagrams has to be based now on an 
auxiliary propagator whose asymptotics coincides with the exact large momentum behaviour and 
is independent of the possibly environment dependent infrared data (masses, condensates) of the 
theory [6]. The exact propagator is expanded around this propagator, which eventually determines 
the counterterms. This analysis becomes unavoidable when also the so-called "basket-ball" diagram 
is added to the set of skeleton diagrams determining the 2PI effective action. It is in Section 6 where 
we extend the analysis to the O(N) model. A central question of these sections will be, how the 
form of sub divergences changes, whether all sub divergences can be expressed as products of a unique 
(possibly environment dependent) function and some explicit cut-off dependent factors, or one has to 
consider several independent "mixed" products. It turns out that the relevant contribution from the 
"basket-ball" diagram can be decomposed into the sum of terms from which some are proportional 
to powers of the background field, some to the finite part of the tadpole integral, therefore no new 
type of subdivergence shows up. This decomposition is determined by an integral kernel solving a 
Bethe-Salpeter type equation, which introduces quite naturally a resummation for certain 4-point 
functions. 

As a consequence, the classification of the divergent pieces remains unchanged relative to the 
classes found already in the 2PI-Hartree approximation [24J: one finds overall divergences propor- 
tional to the zeroth and second power of the background field and sub divergences proportional to the 
finite part of the independent tadpole integrals. Their coefficients are required to vanish separately. 
In general, the number of renormalisation conditions arising from this requirement is larger than 
the number of independent counter couplings. Renormalisability of the 2PI approximation means 
beyond the determination of all resummed counter couplings also the demonstration that these sup- 
plementary equations do not contradict those which are used for extracting the counterterms. 

The summary of our findings is given in Section 7 where we shortly outline possible approaches 
to renormalised (cut-off independent) numerical calculations. Some technical details relevant to the 
discussion presented in the main text are given in two Appendices. 

2 The models and some notations 

The complete static 0(X 2 ) skeleton 2PI functional in the broken phase of the real one-component $ 4 
model, characterised by the vacuum expectation value $ = v is of the following form |13j : 

V[v,G] = imV + - l -Jl\nG-\p) + D-\p)G{p)] + ^ (^Wf 

-^W / 'G(p)G{k)G(p - k) ~ f f fG(q)G(p)G(k)G(p + k + q)+V ct [v,G),(l) 
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where %D 1 (p) = p 2 — m 2 — \v 2 /2 is the tree-level inverse propagator. The counterterm functional 
is given by 



V ct [v, G] = l -5mlv 2 + ^ + l [5ml + f o(p) - f Jp 2 G{p) + 



5Ac 
8 



G{p) 



(2) 



Note, that there is no counterterm diagram involving the "setting-sun" or the "basket-ball" diagrams. 
Such diagrams would be necessary only when going beyond the present level of the 2PI truncation, 
as will become clear in Section 4. 

The generalisation of the effective action to the iV-component scalar field with O(N) symmetry 
in the symmetry breaking background is expressed with the following representation of the tree-level 
and exact propagators: 



\ 2 - m 2 - x 4 

QN 



Xv" 2 

\ r ab ' ab) gjy 



p; b , c ab {k) = gmp:, + gmpi. (3) 



Here P° b 



V Vh P\ 
VaVbi 1 a b 



5ab — v a Vb are the relevant orthogonal projectors. v a , a = 1 . . . N is the 
component of the unit vector giving the orientation of the symmetry breaking constant background. 
In terms of the two propagating modes the 2PI-action has the form: 



V[v a ,G a ,G w ] = ^m 2 v 2 a + 



[v 2 ) 2 

24N y a> 



- I [(N-l)\nG-\k) + \nG-/(k) + (N-l)D;\k)G w (k)+D-\k)G a (k)] 



+ 



A 



24iV 
i\ 2 



36iV 2 a 
i\ 2 



UAN 2 



3^G a (k^j + (N 2 -l)^G 7T (k)^ + 2(N-1) fG^{k) JgM 
3G a (k)G a (p) + (N — l)G n (k)G n (p)] G a (p + k) 
(N - 1) ((jV + l)G n (k)G n (p) + 2G a {k)G a {p)) G„(q)G n (k + p + q) 





k J p J q 



+ZG u {k)G a {p)G a {q)G«{k + p + q) 



+ V ct , 



(4) 



where iD 7r 1 (p) = p 2 — m 2 — \v 2 / (QN),iD a 1 = p 2 — m 2 — \v 2 / (2N) are the tree-level pion and sigma 
inverse propagators. In this case the counterterm functional reads 

Vet = + ^{vlf + \ (5ml + ^) jf (G a (k) + (N — l)G n (k)) + ^ Jc^k) 



+ 



2AN 



5K 



((N - l)G n (k) + G a (k)) 



25\z 



(N-l) 



G n (k) 



+ 



G a {k) 



5Z 



5Z 



~ ir / P G*( P ) ~(N- 1)— / P 2 G W (P)- 



(5) 



The introduction of the two quartic counterterms 5\q and 5Xq corresponds to the two independent 
O(N) -invariant combinations G a bGb a and G aa Gbb- Similarly 5m 2 ,, 5\^, and 5Af are introduced to 
cancel divergences emerging from the two terms, 5 a bG a b and v a VbG a b, appearing in D~ b Gb a - 
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In the one-component model divergent contributions to the equation of state SV/Sv = and to 
the propagator equation 5V/5G{k) = will be expressed with the divergent pieces of the following 
integrals: 

T[G]:= [ G(k) = T div + T F [G], (6) 

Jk 

I(p,G):=-i [ G(k)G(k + p)=I div + I F (p,G), (7) 
Jk 

S(p,G):=-i [ [G(k)G(q)G(k + q + p) = S div (P,G) + S F (p,G). (8) 

Jk Jq 

For the iV-component case one introduces correspondingly T a , I a p, S a /3 7 , where the indices can take 
the values a = {n,a}. Separation of the divergent pieces is explicitly presented in the Appen- 
dices. Without introducing any particular notation, the finite parts of these integrals would include 
temperature and density dependent contributions. 

The non-perturbative construction of the counterterms will be presented in two stages. First, we 
restrict the discussion to the two- loop truncation of the effective actions given in ([1]) and (j4j). In 
this case the asymptotic behaviour of the propagator remains as it is at tree-level. In a second step 
we shall discuss the effect of adding the "basket-ball" diagram, which changes the large momentum 
asymptotics. We believe that the structured presentation makes more clear the algorithm underlying 
the somewhat complex formulae and facilitates its prospective numerical implementations. 



3 Analysis of the two-loop truncation: one component real 
scalar field 

In this truncation one omits the contribution of the last but one term of ([I]). The equation for the 
propagator reads as 

l G-\p) =p 2 - m 2 -Sm 2 2 -^(\ + 5\ 2 )v 2 - ^(A + S\ )T[G} - h 2 v 2 I{p, G). (9) 

The large p asymptotics of the propagator is unchanged compared to the tree- level one, and we 
parametrise the finite self-energy of the exact propagator as 

iG _1 (p) =p 2 - M 2 - n(p). (10) 

The splitting of the finite self-energy into the two pieces 

M 2 = m 2 + K 2 + ±T F [G), n(p) = ^-v 2 I F (p, G), hm _ , (11) 

I I 2, p^oo p L 

separates the momentum-dependent part of the finite self-energy (II (p)), coming from the setting-sun 
integral of the effective-action (CQ) and the momentum-independent contribution (M 2 ) which comes 
from the other two- loop piece of Vfv, G]. Although this latter term formally corresponds to the 
2PI-Hartree truncation discussed in our previous paper [21], its self-consistent determination takes 
into account the effect of II (p). 
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The divergence cancellation condition which determines the counterterms 5m 2 ,5X 2 , and 5X$ is 
obtained by subtracting (TTUT) with the parametrisation (TTTT) from ([9]) and has the form: 



'121 



= 5m\ + l -5\ 2 v 2 + + 5A )T div + ~5\ T F [G] + ~AV/ div . 



The equation of state reads as 



v ( m 2 + 5m 2 + i(A + 5A 4 )t; 2 + J(A + 5A 2 )T[G] + -A 2 ^, G) 
o 2 o 



The cancellation condition for the divergent pieces of this equation amounts to 

= 5m 2 + \5\ 4 v 2 + ~(A + 5A 2 )T div + ^A 2 T F [G] + iA 2 S' div (0, G). 
o 2 2 o 



(13) 



(14) 



From this equation one should determine 5ml, SX^ making use of the previously determined expression 
for 5X 2 . 

The regularised expressions of the counter couplings as determined from (IT21) and (|14p give regu- 
larisation independent solutions when used in the regularised equations @ and ([TBI . In our discussion 
we shall have in mind cut-off regularisation of the divergent integrals. Since the divergent part of 
these integrals is separated using propagators which are not sensitive to the infrared features of the 
theory, their entire zero temperature part will be subtracted. 

With help of Appendix A one finds the following expressions for the divergent pieces of the basic 
integrals: 



T 



drv 



Tf + (M 2 - M' 2 )l^> + ^AVT £ 



,(0) 



2„,2rp(I) 

d ' 



Idiv — T, 
3 



(0) 



5 div (0, G) = S PV (0) + 37fl> + 3(M 2 - Ml) ((Tf ) 2 + if ) + |aV {if if + 1^) , (15) 

The cut-off dependent integrals Spy(O), if, if, if, and if ^ are defined in Appendix A. Mq is 
a free parameter introduced in the auxiliary propagator (j§T|) to ensure the infrared finiteness of all 
integrals. Exploiting the definition of M 2 given in ffTTj) when using (|15|) in (fT2"j) and (|14p . one separates 
the coefficients of v°,v 2 and of Tp[£r] exactly as it was done when we discussed the renormalisation 
of the 2PI-Hartree approximation [24J. The following six equations arise: 



= 5m\ + - (A + 5X ) \lf + (m 2 - M 2 )Tf 



= 5X 2 + ^A(A + 5X ) (if + AT d (/) ) + A 2 T ( 



2 

0= 5A + ^A(A + 5A )Tf , 



2^(0) 
d ' 



= 8ml + i(A + 5A 2 ) [if + (m 2 - M 2 )Tf ] + \x\rn 2 - Af 2 ) ((Tf ) 2 + if) + ^^(0), 
0=5X, + h(X + 5A 2 + A 2 Tf ) (Tf + Alf ) + ^A 3 (rf + Alf 



= 5X 2 + ^A(3A + 5A 2 )Tf + ^A 3 ((Tf ) 2 + Tf 



(16) 




\+6X„ M 2 _ M 2 

CO (") 



Figure 1: Illustration of the appearance of the same radiative structure producing (i) the diver- 
gence denoted with Tjf in the equation for the inverse propagator and (ii) the subdivergence of the 
equation of state. Internal lines represent the auxiliary propagators Gpy. Wiggly line represents the 
background v. Letter 'F' refers to the finite bubble contribution Ipv,F- 

where the first three equations come from (1121) . while the last three from (tHJ) . The expression for 
5Xq determined from the third equation is identical to the expression which we would have obtained 
in 2PI-Hartree approximation. If one determines from the first five equations SttLq, 5m\, <5Ao, 8X2, 
and 5A4, then the sixth serves for a consistency check of the renormalisation. It can be verified by 
direct substitution of the previously obtained counter coupling expressions. The reason behind this 
consistency is that when the exact propagator is expanded around the auxiliary propagator the same 
integral appears both in the momentum independent part of the self-energy and in the contribution of 
the setting-sun to the equation of state. In the former case the divergent part of this integral denoted 
with T { d I] appears as an overall divergence proportional to v 2 , which can be seen by looking at the 
last term of in (TT5l) . In the second case it appears also as a subdivergence, that is proportional 
to T F [G] (see the third term of Sdi v i n (Tl5l) )- This structural correspondence is illustrated in Fig. CD 
We remark that the results in (fl6l) obtained with the present method reproduce results announced 
by Arrizabalaga and Reinosa in a conference poster [TI] (explicit expressions for the counterterms 
were given only on the poster, but not in the Proceedings). The use of these expressions in @ and 
( TT3l) makes these equations renormalised (cut-off independent). 

4 Adding the "basket-ball" diagram 

Keeping all terms in the 0(X 2 ) truncation of the skeleton decomposition of the 2PI potential the 
finite self-energy contains three terms: 

iG~ x {p) = p 2 - E(p), S(p) = M 2 + n (p) + II 2 (p), 
M 2 = m 2 + l -Xv 2 + ±T F {G], n (p) = \x 2 v 2 I F (p), U 2 (p) = jS F (p). (17) 

The cancellation of the divergences of the propagator equation amounts to the equality 

= 5m 2 2 + l -5X 2 v 2 + h\ + 5A )T div + ~5XqT f [G] + ]-X 2 v 2 I div + ]:X 2 S d M ~p 2 SZ. (18) 
2 2 2 2 

The new feature of the analysis is the fact that the large momentum behaviour of the propagator 
changes due to the setting-sun contribution. The new leading asymptotics is accompanied also by an 
infinite wave function renormalisation. The situation is complicated by the fact that the asymptotics 



7 



can be determined only in the course of the numerical solution of the renormalised propagator 
equation (the equation of state has no influence on the asymptotics) . Nevertheless, one proceeds by 
treating separately the leading asymptotics not depending on any infrared parameter of the theory: 

n 2 (p) = n (p) + n 2)0 (p) + M P ), iim ^ ~ (\n P f, iim - o, n r (p) ~ v -\ 

p— *oo p z p^oo p z 

(19) 

The splitting is somewhat arbitrary and this freedom can be exploited to find simple-to-handle ex- 
pressions for the counter-couplings. There is a corresponding splitting of the setting-sun contribution: 

S F (p) = S a>F (p) + 3 , F (p) + S r (p), (20) 

where S a (p) is defined as the setting-sun diagram computed with a conveniently introduced, infrared 
stable auxiliary propagator which asymptotically coincides with the exact propagator: 

iG-\p) =p 2 - Ml - U a (p), G(p) = G a (p) + 5G(p) 7 8G(p) ~ 0{p~ 4 ). (21) 

Weinberg's theorem [23] guarantees that the self-consistent solution of this equation for the self- 
energy (inverse propagator) behaves asymptotically as ~p 2 (lnp) /3 , where j3 can be determined only 
numerically. 

Having introduced all the relevant notations we proceed now with the analysis of the divergences. 
The momentum dependent divergent part of S a (p) is removed by the wave function renormalisation 
counterterm: 

sz = YW SfV{p) - (22) 



(23) 



Again, the renormalisation of the bubble integral is the simplest: 

I(p) = I div + I F (p), I div :=Tj°>, I F (p) = -i [ [G(k)G(k + p)-G 2 a (k)]. 

Jk 

The tadpole integral is analysed with help of the identity given in (11021) 

T[G] = 7f) + (M 2 - M 2 )Tj°) + \\ 2 v 2 T^ - ij G 2 (A;)n 2 , (A;) + jf[G\. (24) 

The finite contribution denoted with Tp\G] consists of two pieces. The first is the contribution of 
G r {p) and the second one is the contribution of I F {k) — I a>F (k), which arises when in U. (k) one 
replaces I F (k) by I a , F (k) to obtain the divergence Ta \ One finds the definitions of T^jT^ and of 
Ta^ in Appendix B. 

The determination of n 2j o(^) is needed to understand the contribution of the last explicit integral 
term in the expression of T[G] given in (1241) . It implies the divergence analysis of S(p). With help 
of the decomposition G{p) = G a (p) + SG(p) (see eq. fl 102ft ) one can disentangle the pieces of the 
setting-sun integral defined in (jBJ), which potentially contain divergences: 

S( P ) = S a (p) + 3 [ 5G(k)I a (k + p) + Si 1 \p), 5«(p)~p- 2 . (25) 



The momentum dependent part of the first term determines the wave function renormalisation, as 
we have already indicated, and also U a (p). The last term gives a contribution to U r (p). 
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One identifies <So,f(p) with the contribution of the integral in ( J25l) which is finite and does not 
vanish asymptotically for large p. Using the expression of 5G given in ( 11021) . simple power counting 
shows that its G r (k) term gives a second contribution to IT r (p). Moreover, in the expression of 5G(k) 
one can replace Ip{k) by I a ^p(k), since the difference gives a third contribution to U r (p). Taking 
in (121)1) the finite part of the asymptotic bubble I a , one observes that the integral has a logarithmic 
divergence which is captured by changing I a ,F(k+p) to I a ,F(k). It is their difference which contributes 
to So,f{p)- To improve the convergence of the integral in the ultraviolet, one exploits the symmetry 
property of the integrand and uses instead of I a ,F(k + p) — I a ,F(k) the kernel 



A 2 r 

K(p, k) = — y a A k + p) + !aA k - p) - ii a A k ) 



(26) 



which is symmetric under the reflection of k — > —k. Note, that the kernel is not symmetric upon the 
interchange of k and p, that is K(p, k) ^ K(k,p). This is to be contrasted with the symmetric kernel 
of Ref. [B], denoted with A(p, k) which, however, is not finite. 

With the above steps the following inhomogeneous integral representation is obtained for Sq^f{p)'- 



X 2 
o 



-- j Gi(k)K(p,k) 



A 



A 2 



M - M 2 + —v'I a Ak) + —SoAk) 
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(27) 



This can be "inverted" iteratively to take the form 
A 2 



6 ' y 1 2 



\(M 2 -M 2 ) (-ijv(p : k)G 



(k) 



+ 7 AV 
4 



-i jT(p,k)G 2 a (k)I a , F (k)y (28) 



where the kernel T(k,p) fulfills the Bethe-Salpeter-like equation 



r{ P ,k) 



k(p, k) - - 



G 2 a (q)K(p,q)T(q,k). 



(29) 



The virtue of the representation in (1281) is that it decomposes So^(p) into the sum of terms 
proportional to v°,v 2 and T F [G}. It is not a solution for S 0: p, since T F [G] contains it implicitly, 
but it means that no new independent finite quantity appears in 112,0 (p) relative to the analysis in 
Section 3. In other words, the sub divergences at this level are found again by looking at the coefficient 
of Tp[G] in the divergence produced by n 2i o(/s) = A 2 S'o i f(p)/6 through the integral in (E 

Using the expression of M 2 from (1T71) in (1281) one obtains 



n 2) o(p) 



1 ' ' M 2 )r (j9) + j [AT (p) + \ 2 Up)} + ^T F [G}T (p), 



[m 



where for the finite integrals we introduced the following shorthand notations: 



r„(p) = -i / T( P ,k)G 2 a (k), 

Jk 



ri(p) 



T{p,k)G 2 a {k)I a Ak). 



(30) 



(31) 



The coefficients T (p) and I\(p) multiplying the quantities sensitive to the infrared parameters (e.g. 



Mq,v ,Tp[G]) are fully determined by G a (p) and have logarithmic large-p asymptotics! 



With the expression of n 2 (p) given in ( 1301) one returns to the expression of the tadpole integral 
241) . Introducing for the divergences of the remaining integral term the following notations: 



Do = -i I G 2 a {k)T Q {k) , D v = -i [ G 2 a (k)r x (k) 
Jk dlv Jk 



di\ 



(32) 



9 



one obtains 



T[G] = + (m 2 - Ml) (t® + ±D ) + \x 2 v 2 [if + AT« + 1(D + AD X ) 



+ 



T F [G]. 



(33) 



Note, that the finite parts of the integrals given in ( l32l) are added to T^fG] of ( l24l) to give Tp[G] 
The analysis of the divergences coming from the tadpole integral is now complete. 

In order to find the expressions of the counter couplings one needs also 5div Using ( 11021) in (1231 
one readily obtains: 

6 



SUP) ~ ^P 2 5Z = S a , div (0) + 3T a (0) (T[G] - T&) 



3* / G 2 a (k)I aiF (k + p) 



M - Mq + -X 2 v 2 I F (k) + n 2l0 (fc) 



div 



SaMV + 3T<?\T[G) - TP) + 3(M 2 - M 2 )Tj J ) 
+ 3^2^(1,2) _ g. f G l (k)IaAk)U2 o(k) } 

I J k div 



(34) 



where T[G] is given in f[3"3"j) and the divergences of the last term are expressible using ([3"U]) . In going 
from the first to the second form of we are allowed to replace both Ip(k) and I a ,F(k + p) by 
Ia,F{k) because the differences give finite contributions. It is remarkable that again the infrared 
quantities (m 2 , Mq , v 2 , Tp[G}) appear linearly with coefficients fully determined by G a - Therefore 
the zero temperature part of the integrals determining these coefficients will be fully included into 
the counterterms! 

Having the expressions of Jdi v , T div , and S^v (see (I2"3"j) . (I3"3l . and (1511 ). one substitutes into the 
last one the definition of M 2 from llTTj) and one proceeds to the decomposition of the divergence 
cancellation condition for the propagator equation into the three separate equations corresponding 
to the pieces proportional to v°, v 2 , and Tp[G], respectively: 



= 5m 2 2 + 1 (A + SX ) [Tf + (m 2 - M*)lf>] + y5 a , div (0) 



+ K - M 2 ) { £ 



+ ^(A + 5A + A 2 ri°))D } 



(35) 



= <5A 2 + -A(A + 6X ) (T a (°) + AT«) + A 2 ^ 
1 



4* 



(T (0) ) 2 + r W +A(T (0) r (/) +T (/, 2)) 

(A + 5X + A 2 TW) (D + ADi) + A 2 (D 1 + AD 2 ) 



= <5A + -A(A + 5Ao)ri ) + A 2 



Ti 0) + ^A((Ti°)) 2 + ry)) 



(A + 5A + A 2 T a (0) ) D + X 2 D 1 



(36) 



(37) 
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The new quantities appearing in these equations are the following: 



£>i = -i / G 2 a (k)I a , F (k)T (k) 



div 



D 2 = -i G 2 a {k)I a , F {k)Y 1 (k) 



div 



(3? 



Note, that D\ ^ D\ because the kernel K(p, k) given in (1261) is asymmetric in its arguments. It should 
be emphasised that at 0(X 2 ) one has 5\q = 5X 2 = —3X 2 Ta°^ /2, as dictated by perturbation theory. 
Also, one can identify in these equations the contribution from each skeleton diagram separately. 
The condition for cancellation of the overall divergences in the equation of state (j!4p determines 



5ml an d ^4 : 



= 8ml + £(A + SX 2 ) [T® + ( m 2 - M 2 )T®] + ^S a>div (0) 



+ (m 2 - M 2 



(T (0) )2 + T y ) + 1^ 



+ \(X + 5X 2 + X 2 tM)D } 



= ^ + ^A(A + 5X 2 ) (Tf) + ATf )) + ^A 3 [(T^) 2 + T« + A (T^ + T^) 
+- A X [ (A + 5A 2 + A 2 Ti°») (D + ADi) + A 2 p x + XD 2 



(39) 



(40) 



In addition, the subdivergence cancellation in the equation of state leads to a consistency check 
between SXq and 5X 2 . The simplest way to obtain it is as follows. Denoting by Tdi v (Tp) and 
So,(Hv(Tf) the contribution from the tadpole and the setting-sun integrals proportional to T F [G], 
the subdivergence cancellation in (fl4l) reads as 



hx + SX 2 )T div (T F ) + ~5X 2 T F [G] + h 2 S 0Aiv (TF) = 0. 
Z l o 



(41) 



One should compare this with the condition for subdivergence cancellation in the propagator equation 
(USD: 

1 1 1 _ 

(42) 



^(A + 6X )T div (T F ) + ^A T F [G] + ^^(TV) = 0. 



Using ( 1331) . from the comparison it follows that 

(<5A - SX 2 ) 



1 + 1 (if + ift) 



0. 



(43) 



It is not obvious, when one looks at (IHUj) and (1371) . that this relation is fulfilled, since when one takes 
from these equations the difference of the two quartic counter couplings, one finds 



5X — 5X 2 



A a 



Ti 0) + ^0)] _1 {if > +\d 1 + X + ^) 



(44) 



Comparing (J43j) and fj44|) one has 



5A = 5A 2 , 



(45) 
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and the requirement of the vanishing of the expression in the curly bracket of (144j) . This condition 
relates several cut-off dependent integrals which can be (and should be) checked only when the 
numerical solutions for G(p) and for the Bethe-Salpeter-like equation (1291) are found. 

With a fully analogous line of thinking one can investigate the consistency of the cancellation 
conditions of the divergences proportional with v and v 2 in the propagator equation and the equation 
of state. Using ( 14*31) the consistency requirement leads to the relations: 

5m 2 2 = 5m 2 , 5X A = 3{5X 2 + A 2 T a (0) ). (46) 

These two relations can be obtained also by comparing (1351) with (|39|) and (|36|) with (f40j) . Since 
5Xq = — 3A 2 Ti°' ) /2, to 0(X 2 ) all three quartic coupling counterterms are equal, as is expected from 
perturbative renormalisation. 



5 Analysis of the 0(N) model at the two- loop level of the 
truncation 

The broken symmetry phase of the O(N) model represents a non-trivial test-ground for checking the 
consistency of our counterterm determination. It has two independent excitations: the (N — l)-fold 
degenerate "pion" multiplet and the more massive "sigma" . Therefore there are four renormalisability 
conditions for each of the two propagators. Since from the counterterms introduced in Section 2 
only five enter the two propagator equations, one has to check if the remaining three conditions 
are fulfilled. Although the algebra of this calculation is quite challenging with help of symbolic 
computer calculations we could check that the extra conditions do not impose any new relation on 
the counterterms, being identically satisfied. 

The equations for the inverse propagators iG~ x {p) = p 2 — Sg(p) and iG' 1 ^) = p 2 — £qG°) ^ eac ^ 
to the following equations for the respective self-energies: 

EJ(p) = m 2 + 5m 2 2 + (A + 5X^ + ± [(N + 1)A + (N - 1)6** + 25A B ] T„ + i (A + 5X$) T a 
X 2 v 2 

MP), (47) 



9A^ 2 

v 2 1 TV — 1 

EJ(p) = m 2 + 6m 2 2 + (3 A + 6X* + 25 X* ) — + — (3 A + 5X$ + 26 X*) T a + (A + 6X$) T n 

\ 2 2 

Here one ought to recall that the pion propagator determined self-consistently in the 2PI approxima- 
tion does not obey Goldstone's theorem e.g. Sq(p = 0) 0- There exists a proposition at Hartree 
level |26j to modify the 2PI functional to cure this deficiency, which one might attempt to generalise 
when including 0(X 2 ) skeleton diagrams. The method of counterterm resummation presented here 
would work equally well for such modified 2PI functionals. 

When looking at (I48p . one can make the following identifications of the respective finite tadpole 
and bubble contributions: 
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Mi »r + J^v 2 + ^[3T a>F + (N-l)n >F ], W (p) 



x_ 

2N~ ' 6N l ~"' r ' v " -'-""J' 18N 2 
The conditions for the renormalisability of the two self-energies and £q(p) read as 



2„,2 



A 2 w 



[9C(p) + (iV-l)C(p)].(49) 







+ Jn V ' + 6iV ^ + 1)A + (iV " 1)5A ° + 25A ^ T ^ div + 6iV (A + 5A ° ) Tct 



div 



5Xi 



2„,2 



A 2 w 



[(JV - l)5\i + 25A B ] 7> + -^T C , F + — T d 



(0) 



6N 

5m\ + (5A^ + 25Af ) 
1 



6AT 



-(A + *A 4 )T» 

,div 



(5A^ + 25\*) T a p + 



AT 



-*A^V + ^(iV + 8)Tf. 



(50) 



(51) 



For the three bubble integrals we used the analogue of fl92l and the fact that their divergence does 
not depend on the kind of fields which constitute them. 

In analogy with the case of the one- component scalar field, one can write down the explicit 
expressions for T^div and T^ div expanding, as in ([53"]) . the exact propagators in the tadpole integrals 
around the common auxiliary propagator. One obtains 



- 7r,div 



Tf + (M 2 - M 2 )lf , ()X2 
^,i I v = Tf + (M CT 2 -M 2 )Tf + 



2„,2 



X 2 v 



__j rpO-ir(I) 

!)• 

X 2 v 2 



2N 2 



\ 2 2 

T^ (7) + — ^(A^-l)T" 
18N 2 



7T7r(7) 



(52) 
(53) 



Although, the divergence involving integrands proportional to the finite bubbles Jfg does not depend 
on the kind of fields which constitute them, in the notation above we explicitly keep track of the 
bubble which produces the divergence. The upper double index gives the two propagators constitut- 
ing the bubble, while lower index refers to the self-energy they contribute to. This distinction will be 
exploited at the end of this section to generalise, to the present case, the diagrammatic interpretation 
of the consistency of the renormalisation conditions given at the end of Section 3. 

Substituting f[5"2l and fl5~3l into the renormalisation conditions, one uses the defining equations 
of M 2 and Ml and requires again separate vanishing of the coefficients of the divergent pieces 
proportional to v ,v 2 ,T ntF , and T atF , respectively. For the pion propagator one finds: 



= 6m 2 2 + — [(N + 2)X + N6X 



26X*} 



T, 



(2) 



[m 



= 6X% + [(N + 1)A + (N - 1)6)4 + 25A^] 



A 



T (o) + 



MM 
A 2 



(54) 



+(A + 5A 



A 



2N d + 



2N 2 



rpo-a(l) 
1 a,d 



18N 2 



QN a ' 9N 2 
(N - 1)T* 



T. 



a,Tr(I) 



7T7r(7) 

d 



2A 2 



+ 3N Td 



(0) 



(55) 



A 



= {N- 1)6 X^ + 26 X* + — {(N-1)[(N + 2) (A + 6X$) + 2(A + 6X*)} + 4(A + 6X*)} if ,(56) 







6Xt 



X 



— [(N + 2)(X + 6X£)+2(X + 6X*)}T c 



(0) 



(57) 



The last two equations form a closed set for 6Xq and 6X$ . Their expressions can be used in the first 

"2 u/y 2 ■ 



two equations to obtain 5m\ and 8X 2 
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In case of the sigma propagator, the corresponding set of renormalisation conditions reproduces 
exactly the same equation for the mass-renormalisation 8m\. The condition for vanishing of the 
divergent coefficients of T n< p and T a p produces formally different expressions when compared to 
those in f[5*6"]) and (|5"T|) . but their solution for <5Aq and 5Xq coincides with the solution of f[5*B"]) and 
( 15"T|) . The only new condition arises from the vanishing of the coefficient of v 2 , which reads as 

= 5X* + 25Af + ^T] 0) [2(N + 8) A + (N + 2) (A + 5X$) + 6(A + 5X$)] 



A 2 



18N 2 



2(X + 5X*) 



(N - l)T™ d (I) + 9T. CT<T{7) 



+(A 



(N-l)(TZ {I) + 2TZ {I) ) + 9TZ {I) 



(5? 



Using the expression of 5X^ obtained from (1531) the equation above determines 5Af . Note, that the 
expressions for 5Xq, 5Xq and bm\ are identical with those obtained in the 2PI-Hartree approxima- 
tion in [21]. The extra term in the effective action, compared to the 2PI-Hartree approximation is 
responsible for the fact that 5X^ ^ SXq and <5Af ^ 5Xq. 

In the broken symmetry phase the equation of state is given by 







m ' + 6m2 ° + 6^ (A + 6Xi)v2 + ^AT ( A + 6X ^ T ^ + ^N ( 3A + 6X2 + 2X ^ Tu 



+ 



A 2 



18N 2 



3^ CTCT (0) + (iV-l)5 CT _(0) 



(59) 



The renormalisation conditions emerging from this equation naturally separate into two groups. The 
vanishing of the divergent coefficients proportional to v and v 2 provides the expressions of 5m 2 , and 
6X4. The consistent renormalisability of the present truncation of the 2PI-approximation depends on 
the vanishing of the sub divergences proportional to T^p and T as p. There are no extra free parameters 
to ensure this. 

In order to explicitly write the renormalisation conditions, one needs the divergences of the two 
setting-sun integrals. To obtain them, one does the same steps as in ( |93l) and derives expressions 
similar to f l94l) for SG a (p) = G a (p) — Gpy{p) and SG v {p) = G v {p) — Gpy(p), namely: 



6G a (p) 
SG^ip) 



iG Pv(p) 



iG 2 PV {p) 



Ml 



Mi 



9C(P) + (N - l)l£(p) 



M, - Mi 



X 2 v 2 
18N 2 



+ G r a (p), 



(60) 
(61) 



The regular parts G r a (p) and G r ^{p) behave asymptotically as p~ 6 . Using G a (p) = Gpy{p) + 5G a (p) 
with a = {a, it}, one finds the following expressions for the two types of setting-sun integrals in 
which the parts contributing to the divergences are well separated: 



Saaa(0) = S PV (0) + 3 / 8G a {k)I p ^{k) + S* f) , 

Jk 

Sow(0) = S PV (0) + / 6G a (k)I PV (k) + 2 / 6G n (k)I PV (k) + S%$. 

Jk Jk 



(62) 



For the bubble I P v{k) we indicated the original two propagators which were replaced by auxiliary 
ones. 
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When substituting ( J60l) and (IBTl) into (l62l) one replaces I^pip) with Ipy^ F {p) (« £ {7r,cr}), since 
the difference gives a finite contribution. The same steps as performed on ( l9Tj) give for the divergent 
parts of the setting-sun diagrams: 



s£s;(o) = s Py (o) + 3Tf ( r CT , F - rr ) + 3(m ( 



,(2) 



Spy(O) 
AV 



vj 0) (/ -/ 



3Tf ) + (M 2 - M 2 )^ + 2(M 2 - M 2 )T, 



M o 2 ) r M (/) + ^(^ + 8 ) T f 2) 

■7T7r(/) 



-lO-7r(7) 



+ 



18iV 2 



(iV+ 12)T d (/,2) . 



(63) 



(64) 



The divergent expressions , T^^\ and T^ l) are the same as those occurring in (1521) and 



<77r(I) 



For the divergent part of the integral containing two bubbles we used the definition (1991) without 
distinguishing which bubble produces it, because this is an overall divergence and does not show up 
in any consistency condition. 

The condition for cancellation of the divergences in the equation of state (151?]) reads as 



= 5m 2 . + — -V + 



5X 4 
6iV 



6N 



1 A N 

-(A + 5X 2 )T 7ridiv + 



6iV 



+(3A + 5A 



6N 



+ (5X$ + 25Xt 



+ 



A 2 



6iV 18iV 2 



3^-(0) + (N - 1)^(0) • (65) 



The renormalisation conditions can be read from this equation using the expressions given in 
( l53l) . ( |63l) . and ( !64l) with the corresponding formulae for M 2 and M 2 given in ( 1491) . The counterterms 
5m 2 , and 5A4 are determined from the vanishing of the coefficients of v and t> 2 , respectively. The 
conditions for the cancellation of the coefficients of T a% p and T nt p are given by 



= ^(^ + 2^ 

A^ 
+ AT3 



AT 



(0) 



36iV 2 



3A(iV + 8) + (N + 2)5X$ + 6SX 



+ 



A 3 (Tj 0) ) 2 
108iV 3 



(5iV + 22) 



-T 

4 °4 



+ 



N 



36 



54 



-T. 



o-7r(7) 



(66) 



<5A^ 



+ 



AT, 



(0) 



6iV 

A 3 
+ 6iV 3 



36iV 2 

}_rpva(I) 

2 CT,d 



A(iV + 8) + (N + 2)5X$ + 25X 



+ 



A 3 (Tf) 2 
108iV 3 



(3iV + 10) 



+ 



N 



-T 



7T7r(J) 



+ 



iV+ 1 

IT - 



T 



o-7r(7) 



7T,«f 



(67) 



Using the solution for 5X 2 and £Af , as obtained from (1551) and ( 1581) . and the solution for SXq and 
5A^ obtained from (156!) and (|57|) . one can easily verify that the above two consistency conditions are 
fulfilled. Note, that the coefficients of T^ , T^ T \ and T^ 7 ^ separately vanish when the solution 
for the coupling counterterms are used. The reason is the same as the one given at the end of 
Section 3, namely the appearance of the same structure in the renormalisation conditions of both 
the propagator equations and of the equation of state. The diagrammatic illustration given in Fig. [2] 
parallels that of Fig. CD Diagrams (a), (b) and diagram (c) correspond to the last two terms of ( 1531) 
and the last term of (1521) . respectively, while diagram (d) and diagrams (e) and (f) correspond to the 
last but one term in (!6"3"j) and the last two terms of the first line of (I6"4"I) . respectively. 
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3\+8\£+28\* 3A+6\£+26\$ \+8\£ A/ 2 -A/ 2 M*-M 2 A/ 2 -J\/ 2 

(a) (b) (c) (d) (e) (f) 



Figure 2: Illustration of the appearance of the same radiative structures producing the divergences 
denoted by T^ z \ T^ T \ and T™} 1 ^ in the equation for the sigma propagator and the sub divergences 
of the equation of state. Solid (dashed) lines indicate that the auxiliary propagator Gpy comes from 
a (pseudo) scalar propagator. For other notations we refer to Fig. [TJ 



6 Analysis of the O(N) model with "basket-ball" diagrams 
included 

At the 0(X 2 ) truncation level of the effective action the finite self-energies are decomposed as 

EM = Ml + nj(p) + nj(p), nj(p) = [s^( P ) + (n + i)^»] , (68) 
K(p) = Ml + w ( P ) + nj(p), nj(p) = ^ B» + (jv - , (69) 

where M%, M%, Hq(p) and Hq{p) are defined in (H9il . Since the infrared sector has no influence on 
the asymptotic behaviour of the self-energies, one can separate a common leading asymptotics in the 
part of the self-energies containing the contribution of the setting-sun integrals (a E {ir, a}): 

n«( P ) = u a (p) + n",(p) + n?(p), iim „ (lnp) i im 51^1 _> 0j n r » ~ P - 2 . 

(70) 

The corresponding decomposition of the finite contribution of the setting-sun integrals reads as 

Sa/Spip) = SaAP) + S awiP) + S£/9/j(p), ( 71 ) 

where a, (3 G {7r, er}. 

Compared to ( 150]) and ( 15~TT) the renormalisation conditions for the equations of G~ l {p) and G" 1 ^) 
are now completed, respectively, by the following terms: 

-p 2 5Z + ^{ v ), -p 2 5Z + ^(p), (72) 

where 



(p) = ^ [5^(p) + (iV + 1)^(P)] , Ef (p) = W [3S^(P) + (iV - 1)S^ 

(73) 

The renormalisation condition (161)1) for the equation of state remains unchanged. Evidently, in (|50|) 
and (15T|) one has to replace by Ti°\ 
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With the same steps which lead from f llOip to (11021) one can write the following representation 
for the propagators 

G a (p) = G a (p) + 8G a (p), 6G a (p) := -iG 2 a (p) [M 2 a - M 2 + n°(p) + Tl^p)] + G«(p), (74) 

where G> (p) behaves asymptotically as p~ 6 . 

Using (1741 . the divergences of the tadpole integrals can be easily separated 



T a = + {Ml - M 2 )Ti°) + ^L{N + 8)T« - i J Gl(k)Ul 
% = + (Ml - M 2 )T(°) + - i jf Gimi 



,(*) 



div 



,(*) 



(75) 



The same procedure which led to (l27"|) provides now an integral equation for the finite piece of 
the setting-sun integrals S^ /3 (p) which grows no faster than logarithmically for asymptotic values 
of p. Combining the corresponding contributions of different setting-sun integrals one obtains the 
following set of coupled integral representations: 



n£, (p) 
nj,o(p) 



18iV 2 J k 
i 



18N 2 



G 2 a (k)K(p, k) [F a (k) + A aa W 2fi (k) + A an Ul (k)} 
G 2 a (k)K(p, k) [F„(k) + A^W^k) + A^W 2fi {k)] 



(76) 



K(p, k) was defined in (126"]) . Here we introduced the matrix A with elements A aa = N + 8, A a7T 
2(N — 1), A na = 2, A nn = 3N + 4, and also the two inhomogeneous terms of the above relations: 



2„,2 



X 2 v 



F a (k) = A aa (Ml - Mq) + A CT7r (M 2 - Mq) + —B^fc), 

AV 



(77) 



= K*{M* - M 1 ^ + A^{Ml - Mi 

with E CT = (N + 8) 2 + 4(iV - 1) and B n = 8(N + 3). 
The solution of ( 1T61) takes the form 



18N 2 



BJ aiF {k) } 



n?o(p) 



18iV 2 



G 2 (fc)r a/3 (p,A;) 



^ 7 (M 2 - M 2 m 7 ) + —^/^(fc) 



(7E 



where 7 G {7r,a}, tt„ 
equation 



1. T a p(p,k) satisfies the following Bethe-Salpeter-like matrix 



(79) 



r>(p, k) = K(p, k)5 af3 - / K(p, q)A ai Y lP (q } k) 



One observes that T aa (p, k) couples only to T na (p, k) and T n7r (p, k) couples to T an (p, k). Exploiting 
that A a . a . — A ncT = A wn —A an , one can verify that the two combinations Qi(p, k) = T aa (p, k)—T na (p, k) 
and Q2(p, k) = T nn (p, k) — T an (p, k) satisfy the same equation 



Qi(p, k) = K(p, k)-c 



187V 2 
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K(p,q)Qi(q,k), 



(80) 



where c = N + 6. This means that Qi(p, k) = Q2(p, k), and consequently there is a relation between 
the elements of the T(p, k) matrix: 

T aa (p, k) - r na (p, k) = T n7T (p, k) - T a7T (p, k). (81) 

Another relation which can be derived reads as 

T aa (p, k) + (N- i)r w (p, k) = r m (p, k) + {N- ly'r^ip, k). (82) 

This is because the combinations on the left- and right-hand sides of the equation both satisfy an 
equation similar to ( IHUl) . but with c replaced by c = 3(A^ + 2). Making use of ( IHTi) in ( 1821) one obtains 

T CTW (p,k) = (N-l)T W(T (p,k). (83) 

We shall see, that the relations (IHTI) and (1H31 play an important role in checking some of the consis- 
tency relations of our renormalisation procedure. 

Having obtained the integral representation ( 1781 for TE^ (p) and ITJ (p), one uses them in ( 1751) to 
write the divergent part of the tadpole integrals with help of the new divergent matrices D^ 1 to be 
defined below: 

T Miv = Ti 2 V + {T^5 ai + D^A M ) (M 2 - M 2 u 7 ) + ^ {r^A aa + DfB p ) . (84) 

The same new divergent quantities appear in the divergent part of the self-energies containing the 
contributions of the setting-sun integrals. One obtains with the procedure which gave ( 1341 the 
following expression: 



Ef (p) = (N + 2)^S B (pK + (^) 2 [(Tfr« + Ti^)) ^ + A aP (tWd* + p 



(T (o ))2 + T (/)U + ( T f)^ + z)M aJ (M 2 - M\) + T^T,A(85) 



A 2 

+ T^^ 

The divergent quantities _D 4 a/3 with z e {0, 1, 2} and Z?™^ are obvious generalisations of those appearing 
in (l32l) and (l38l) with Tq 13 and T"^ given by 

rf(p) = -t^ / r^(p, fc)G 2 (A;), rf (p) = --Lj /" I>(p, fc)G 2 (A;) 

These matrices inherit from T a p some relations between its elements which are analogous to ( JSTl) and 
(J8J): 

D ™ - = n™ - d™, b\ a - b\ u = b™ - d™, (86) 

Df* = (N - l)D™, = (N — \)b\ a . (87) 

In order to determine the counterterms one has to use (J84l) and (l85l) in the renormalisation 
conditions for the two propagators (1501) . (15T1) . supplemented with the terms given in ( 1721 . and in the 
renormalisation condition of the equation of state ( 1651) . In what follows we review, without giving 
the lengthy explicit expressions, how one can obtain the counterterms and the consistency conditions 
similar to 
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The vanishing of the coefficients of v°, v 2 , T a> p, and T„ t p in the two renormalisation conditions for 
G n (p), G a {p) results in eight conditions for the determination of five counterterms: drrt 2 ,, 5Xq, 5Xq, 
SX2, and SX2 ■ Out of them 5Xf and <5Af are obtained from the vanishing of the coefficients of v 2 . 
The vanishing of the divergent coefficients of T a p and T n p in any one of the propagators determines 
SXq and SXq. The condition to obtain the same expressions from the renormalisation conditions of 
the other propagator reduces to the equation 

A 2 (iV + 2) [CiTi 0) - C Q T<p - 3(JV + 2) (C D™ - C X D™ + D^D™ - D™D™)] 

+6N(XCi + NC ) = 0, (88) 

where Cj = (N — l)DJ a — . Using (1871) this consistency condition is identically satisfied. The 
same way, in order to obtain the same expression for 5m 2 , from the vanishing of the coefficients of v 
in either one of the two renormalisation conditions for the propagators one needs 

[X 2 (N + 6)T<® + 6N (A + 5X$)] {D™ + D™ - D™ - D™) 

+X 2 (N + 6) (pi* + D™ - b\ a - = 0. (89) 

By ( |86l) . this is again identically satisfied. As a partial summary we emphasise that the eight condi- 
tions renormalising the propagator equations of the O(N) model in the broken phase are algebraically 
satisfied with the five counterterms, independent of the form of the asymptotic propagator G a . 

Turning now to the equation of state, one can subtract its renormalisation condition from the 
renormalisation condition of G a , and require the cancellation of the term proportional to T at p,T nt p, 
v° and v 2 . We can easily find the following relations: 

5X^ = 5X1 5Af = 5A* 5ml = 5ml 5X 4 = 5X$ + 25X* + ^(N + 8)T^ . (90) 

The first two conditions above impose two consistency conditions on the coupling counterterms if 
one uses their expression as obtained from the renormalisation conditions for the propagators. These 
are complicated expressions, similar to that in (|44p . which now depend also on and relate several 
cut-off dependent integrals. As in the one-component case, these relations can be checked only 
numerically. 

7 Conclusions 

In this paper we have extended the renormalisation procedure developed in Ref. [21] to the case of 
momentum-dependent self-energy appearing when one goes beyond the Hartree-Fock level of trun- 
cation of the 2PI effective action. Two types of momentum-dependent truncations of scalar <3> 4 type 
models were investigated. In the first case the asymptotic behaviour of the exact propagator remains 
unchanged relative to the tree-level, while in the second case the asymptotic behaviour of the exact 
propagator changes. 

Working in the broken symmetry phase we introduced infrared safe auxiliary propagators having 
the right asymptotic behaviour in order to explore the divergence structure of the model. Beyond 
the overall divergences proportional to the zeroth and second power of the background field we found 
only such sub divergences which were proportional to the finite part of the tadpole integrals. This 
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is the reason that no counterterms are needed beyond those appearing at the Hartree truncation of 
the 2PI effective potential. Requiring the vanishing of the coefficients of independent divergences, 
we gave explicit expressions for the counter couplings of the model in terms of some integrals of the 
auxiliary propagator. 

For the first type of truncation investigated, our method reproduces the expressions of the counter 
couplings which were previously reported in the literature for the case of a real selfinteracting scalar 
field, but with no reference to any Bethe-Salpeter equation. Although, the divergence structure was 
known in the literature for the second type of truncation, the explicit expressions for the counter 
couplings, which are needed to make the equations of the theory finite, appear to our best knowledge 
for the first time in this paper. These expressions depend in a non-trivial way on the solution of a 
Bethe-Salpeter-like equation and can be determined only numerically. This was beyond the scope of 
our present study in which we concentrated mainly on the structural aspects of the renormalisation. 

The proposed method treats the conditions for overall and subdivergence cancellation on equal 
footing in the propagator equation(s) and the equation of state. Therefore, the number of conditions 
is larger than the actual number of counterterms. As shown in this paper, these conditions are 
algebraically fulfilled for the first type of truncation, both in the case of a one-component and a 
multi-component scalar field (O(N) model). For the second type of truncation at least some part 
of the consistency conditions can be verified only numerically. It is quite impressive, though, that 
in this case the consistency of the propagator renormalisation can be proved algebraically. In the 
explicit demonstration the internal structure of the Bethe-Salpeter-like matrix equation plays an 
essential role. It is the compatibility of the renormalised propagator equations with the renormalised 
equation of state which leads to those consistency equations which can be verified only numerically. 

The auxiliary propagators applied for the separation of divergences introduce also the renormali- 
sation scale Mq. The explicit dependence of the counter couplings on this scale allows to establish the 
renormalised trajectory of the 2PI approximation in a higher-dimensional coupling space correspond- 
ing to the larger number of independent bare couplings necessary for its consistent renormalisation 
than in the perturbative case. Even when one restricts the beta-functions to leading (e.g. 0(X 2 )) 
order one finds the perturbative beta-function only for the coefficient of the v A term in the 2PI- 
functional. It is interesting to note that beyond the perturbatively higher order contributions to the 
counter couplings also the appearance of the asymptotic propagator G a (p) in the relevant integrals 
produces deviations from the perturbative behaviour of the beta-functions determined in the 2PI 
approximation. 

When passing to the numerical investigations the spirit of our paper suggests to make use of 
the explicit expressions for the cut-off dependent counterterms in the cut-off regulated equations 
(for instance (Q and ffl3l) ) and verify numerically the cut-off independence of their solution. Alter- 
natively, one can prescribe values for a number of physical quantities in terms of the renormalised 
one- and two-point functions. One attempts then the determination of the asymptotic behaviour of 
the propagator from a finite twice-subtracted dispersive representation satisfying the renormalisation 
prescriptions. Before proceeding to the solution of the self-consistent equations in the non-asymptotic 
regime one has to check the consistency relations between asymptotic integrals. Numerical investi- 
gation of these questions essential for the validity of the renormalisation methods proposed for the 
2PI approximation, will be reported elsewhere. 
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A Appendix A: Divergence structures for the two-loop 
truncation of the 2PI approximation 

The divergences of the integrals defined in (jSJ), (J7J) and (jSJ) can be most conveniently explored us- 
ing Pauli-Villars type terms obtained by adding and subtracting in the integrands the auxiliary 
propagator 

iG P l v { P )= P 2 -Ml (91) 

The simplest is to separate the logarithmic divergent part of the bubble integral from its finite 
part. They read as 



'div 



/ G 2 PV {k) =: Tj 0) , I F (p) = -i f [G(k)G(k + p)-G 2 PV (k)] J (92) 
Jk dlv Jk 



where the divergent piece is clearly independent of the external momentum. 
For the propagator appearing in ( 1101) one uses the identity 



G(p) = G PV (p) - iG(p)G PV (p) (V - M 2 + y^Mp)) . (93) 

Solving (I9"3"j) for G(p) and expanding it around Gpy(p) this equality can be rearranged also as 

G( P ) = G PV (p) + 5G(p), 5G(p) = -iG 2 PV {p) (V - M 2 + ^M?)) + G r (p), (94) 

where G r (p) is the regular part of the propagator, which contains at least three powers of Gpy(p)- 
It behaves asymptotically as ~ p~ % and produces a convergent contribution when integrated over 
momentum. 

With help of the representation ( 1941) the divergences and the finite part of the tadpole integral are 
easily separable. After replacing Ip{p) by Ipv^ip) m (jBD (since the contribution of the difference 
to the tadpole is finite) the following partitioning is found: 

T[G] = T[G PV ] + (M 2 - M 2 )Tf - \\ 2 vH JgI v {p)I pv>f {p) + T F [G\. (95) 

Here one introduces the notations 

T[Gp V ] = [ G PV (p) =: Tf\ Tf := -i [ G 2 PV (p)I PV Ap) , (96) 
Jp dlv Jp dlv 

and obtains the first entry in (|15|) . 

The setting-sun integral is rewritten after the replacement G(p) = Gpy{p) + 5G(p) is made for 
each of the three propagators as 



S(p = 0, G) = S PV (0) + 3 / 5G(k)I PV (k) + S { F \ 

Jk 



(97) 



where Spy(0) = — i f k f Gpy{k)Gpy(q)Gpy{k + q) .In the integral above one decomposes Ipv{k) 

1 div 

into a divergent and a finite part. For the piece proportional to the divergent part Ipv,aiv = one 
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uses 5G(k) = G(p) — Gpv{k), while when taking the finite part of the bubble one substitutes 5G(k) 
from ( 1941) . With help of (1961) the following representation for the divergent part is obtained: 



S div (0,G) = S PV (0)+3(T[G]-TP)T^+3(M 2 -M^TP~h 2 v 2 t [ G 2 PV (h)I F (k)I PV , F (k) . (98) 

1 Jk dlv 

The divergent part of the integral above can be obtained by replacing by Ipv,F(k) (the difference 
gives a finite contribution). Then, introducing another cut-off dependent integral 



T. 



one obtains the last entry in ( fl5l) . 



(*.2) 



-i [ G 2 PV {k)I 2 PV Ak) 
Jk 



div 



(99) 



B Divergence structures for the truncation of the 2PI ap- 
proximation at 0(X 2 ) skeleton diagram level 

The analysis goes quite in parallel with that given in Appendix A. One defines the basic divergent 
integrals to be used in the exploration of the divergence structure: 



T (2) 



Ga( P ) 



div 



Ti 0) := ~i ! G 2 M 
J v 



div 



:= -i [ GlWaAv) 



div 



div 



T^ 2) := ~i [ GIWIAP) 
Making use of (fT7|) . ffTUl) and (j2*Tl) . the identity corresponding to flD^j) reads 

G( P ) = G a {p) - iG(p)G a (p) [M 2 - M 2 + n (p) + n 2>0 (p) + n r ( P )] 

Solving (110 If) for G(p) and expanding around G a {p) one obtains 
G'(p)=G a (p) + 5G(p), *G(p) := -zG 2 a (p) 



(100) 



(101) 



1 



M 2 -M 2 + -M F (ri + n 2 , (p) 



+a(p), (102) 



where G r (p) contains also the contribution of H r (p). Power counting shows that for asymptotically 
large values of p one has G r (p) ~ p~ 6 . 

Eq. f)102p is used to separate the divergence of the tadpole integral (see also to analyse the 

divergence structure of the momentum-dependent setting-sun integral and, eventually to obtain the 
integral representation of So,f(p) (e.g. n 2j o(p))- 
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